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We communicate results on correlation functions for tfie spin- 1/2 Heisenberg chain in two particularly impor- 
tant cases: (a) for the infinite chain at arbitrary finite temperature T, and (b) for finite chains of arbitrary length 
L in the ground state. In both cases we present explicit formulas expressing the short-range correlators in a 
range of up to seven lattice sites in terms of a single function co encoding the dependence of the correlators on T 
(L). These formulas allow us to obtain accurate numerical values for the correlators and derived quantities like 
the entanglement entropy. By calculating the low T (large L) asymptotics of u we show that the asymptotics of 
the static correlation functions at any finite distance are (1 /L^) terms. We obtain exact and explicit formulas 
for the coefficients of the leading order terms for up to eight lattice sites. 



Until about ten years ago it was widely believed that it 
would be practically impossible to calculate the lattice cor- 
relation functions of the Heisenberg chain explicitly, despite 
the fact that its Hamiltonian 

L 

■Hat = J^(7jCTj+i (1) 
i=i 

is one of the best known among those integrable lattice sys- 
tems whose spectrum can be exactly calculated by means of 
the Bethe ansatz 1 1 1. Starting from about the year 2000, how- 
ever, our understanding of the model has changed profoundly. 

The recent progress is due to several different but related 
methods, like the representation theory of quantum algebras, 
the algebraic Bethe ansatz, or the advent of new functional 
equations, which first of all led to the derivation of multiple 
integral representations f2\ for the elements of the density ma- 
trix. In a series of works |3 , 4| it was then shown that the cor- 
relation functions of an inhomogeneous version of the model 
can be represented by algebraic expressions of a single func- 
tion of two variables. 

An exponential formula for the reduced density matrix was 
obtained in (5). This formula will be of particular importance 
below, since it is also valid for finite temperature or finite 
length as was conjectured in 161. In Q a fermionic structure 
on the space of local operators was identified from which a 
generating function of all local correlation functions was ob- 
tained 18| for a very general inhomogenous vertex model in- 
cluding the finite temperature and the finite length Heisenberg 
chain as special cases. Even in this most general situation 
the inhomogeneous correlation functions depend on only two 
functions [^,'9'|. Apart from these important theoretical find- 
ings many results for the ground state correlation functions 
in the thermodynamic limit, but also a number of finite tem- 
perature or finite length correlation functions were obtained 
explicitly, see e.g. 1 10 1 and the references therein. 

The aim of this paper is twofold. First, we present re- 
cent exact results for the temperature and size dependence 
of two-point correlation functions. Of particular interest are 
system parameters T and L in the regime of conformal field 
theory (CFT) and their influence on short-ranged correlators 



that are -strictly speaking- outside the domain of CFT. Fur- 
thermore, we present results for the entanglement entropy of 
sub-segments of an infinitely long chain in the entire temper- 
ature window from to oo. Second, we intend to explain the 
necessary computations in the light of a recent understanding 
of the subject 1 11 1 based on 'discrete functional equations' 
which we believe is physically most natural. 

All information about the coiTelation functions of operators 
acting on n successive sites of the chain are encoded in the 
reduced density matrix Z?„. The matrix elements of Dn are 
denoted by i'^i ' ct^' correspond to the expectation val- 
ues of operators |(Ti, (T„) (/ii, /i„|. 

A quantum mechanical system at finite temperature T may 
be viewed as a classical statistical system on a cylinder of cir- 
cumference (3 = 1/T. In this way, the density matrix _D„ on 
n successive sites for temperature T is obtained in a suitable 
limit of the six-vertex model density matrix on a rectangular 
lattice of unbounded width and finite height . Each row cor- 
responds to an imaginary time slice of height r = /3/N, the 
continuous time limit is obtained in the Trotter limit iV — > cx). 
For our purposes it is convenient to have independent heights 
Ti, ...,Tjv (resp. spectral parameters ei,... placed on the 
horizontal fines) under suitable conditions like Tj — 0{1/N) 
and J2j = P- 

We introduce independent spectral parameters xi, a;„ on 
the veitical lines coiTesponding to the n sites picked for the 
definition of The density matrix of the generalized prob- 
lem now depends on the Xj \ and is denoted by _D„(xi, a;„) 
with matrix elements D^j'" -'^^^ (xi, a;„). The full func- 
tional dependence will be solved, the subsequent specializa- 
tion of the arguments yields the physically interesting data. 
The reason behind the success of this solution strategy is the 
analyticity of the object as a function of the spectral param- 
eters X. This is a consequence of the integrabilty of the sys- 
tem: the eigenstates of the column-to-column transfer matri- 
ces T{x) do not depend on x and hence the eigenvalues and 
other objects as functions of x do not show the otherwise un- 
avoidable root singularities. 

As an example of this solution strategy we like to note re- 
sults for the leading eigenvalue h.{x) of the transfer matrix 



2 



T{x) which will be of use below 



logA(a;) = eo(a;) + 



log(l + a(2/)) 



2t: Jc [x-y- i){x - y) 



dy. (2) 



Here, C is a narrow closed contour around the real axis, eo is 
some elementary function containing the -parameters. The 
auxiliary function a{x) satisfies the non-linear integral equa- 
tion 

loga(.) = ao(.)-i/"'$^±^d,, (3) 
TT Jc l + {x- yY 

where ao(a;) is an elementary function of x and the Cj, which 
for finite temperature, in the Trotter limit N ^ oo, takes the 
formao(a;) = 2J/3/(x(x + i)). For the ground-state and finite 
size it is simply aQ{x) = L log fq^jyl- 

The program of calculating the functional dependence of 
Dn{xi, is well understood in the case of zero temper- 
ature. For T = (i.e. (5 — oo) the object Dn{xi, a;„) sat- 
isfies non-trivial linear functional equations in all arguments. 
The most important equation, the 'rqKZ' -equation reads for 
instance 

Dn{xi,...,Xn-l,Xn-l) = A{xi, Xn)Dnixi, Xn) (A) 

where A is a linear operator |j4i| acting in the space of den- 
sity matrices and the Xj may take any value from the complex 
plane. The action of A on Z?„ is of the following type 



(5) 



The functional equations can be solved for the analytical 
function D. The uniqueness of the solution for T = is guar- 
anteed by the asymptotic behaviour for xj — > oo. The find- 
ing of 14J is striking, a two-point nearest neighbour function 
uj{xi,X2) ■— 6 tr D2{xi, X2)SfS2 and a set of 'structure con- 
stants' fn,i,j determine the density matrix £>„ for arbitrary n 

[»/2] / m \ 

E E ( n^(^^p'^'^p) ) /"J,j(^i'2;2, •■•,a:^n) (6) 

m.=0 /,,/ \p=l / 

The summation labels / and J are m-tuples of integers such 
that / n J = and /i < • • • < /„, 1 < /p < < n. The 
structure coefficients fn,i.j{xi,X2, ---jXri) are matrices with 
elementary rational functions of the arguments xi, a;2, a;„ 
as entries. For details see |4J. 

At zero temperature the computation of the correlation 
functions employed two important facts: (i) the functional 
equation Q is satisfied for continuous arguments, and (ii) 
Dn{xi, ...,Xn) depends on the arguments only via the dif- 
ferences Xj — Xi. This is fundamentally different for finite 
temperature, i.e. finite Trotter number N. The main problem 
is that (|4]) no longer holds: on the r.h.s. of (j4]) untreatable 
'correction terms' appear, a serious obstacle so far for treating 
T > 0. 



At this point, the computation of the density matrix for fi- 
nite Trotter number TV on the basis of functional equations 
looks unfeasable. However, for the above six-vertex model 
with N many rows carrying spectral parameters ei, e^v we 
find 1 11 1 a discrete version of functional equations. In detail: 

(I) Equation Q 'rqKZ' holds literally (!) if we restrict a;„ to 
the set {ei, eat}. 

(II) A reduction takes place in the limit of large spectral pa- 
rameter a;„, i.e. D^ixi, ...,x„) D^-iixi, ...,x„_i) (g) Di 
where Di is a single site density matrix of a paramagnetic 
spin. For zero magnetic field this is Di — ^id. 

(III) Analyticity properties: 



A{xi) ■ ... • A{Xn) 



(7) 



where P/^^'---'^^ (xi, x„) is an n-variate polynomial of de- 
gree N in the variables xi,...,Xn and A(x) is the largest 
eigenvalue of the transfer matrix T{x) obtained from (j2|. 

The computational problem consists in finding the polyno- 
mials P^'^i' ■■■^^'(a;i, ...,x„). We have proved |11| that the 
above equations have a unique solution. Hence any expres- 
sion (or ansatz) for Dn{xi, a;„) exhibiting the above listed 
properties realizes the solution. 

We find that (I) is satisfied by (|6]l like in [4J with the same 
structure coefficients fn.i.j provided that u!{xi,X2) is a sym- 
metric function satisfying the 'discrete functional equation' 
(see 1(9] for the most general case) 



Uj{xi,X2) Uj{xi,X2 - 1) 



X2 - 1 



x{x + 2) 2{x^ - l)x{x + 2) 



(8) 



where x := xi—X2, xi is arbitrary and X2 € {ei, eat}. The 
asymptotics (II) is satisfied if Ld{xi,X2) — for xi — ?► oo or 
X2 oo. Condition (III) is satisfied if w is a polynomial (of 
degree N~ 1) divided by A(a;i)A(a;2). These conditions for w 
characterize a unique function. It is relatively straightforward 
to see that the conditions for uj are satisfied by the following 
expressions 



a;(a;i,a:2) := ^ + ^{{xi - X2) 



■>p{xi,X2) 



1 



dy 



2 l)'ip{ixi,ix2) 
G{y,xi) 



Jc ^ + a-iy) (y - X2)iy - X2 - i) 
where the function G satisfies the linear integral equation 



(9) 



G{x,xi) = 

dy/ir G{y,xi) 
{x - xi)ix - xi-i) ' 1 + 1 + (a; - 2/) 



1 



(10) 



and the function a{x) was introduced in ([s]). 

Note that the expression in (j6]l contains rational functions 
as prefactors. The poles of these coefficients are cancelled by 
zeros appearing due to the pairwise cancellation of terms in 
the sum. Consequently, the only poles on the r.h.s. of (|6]l are 
those occurring in the product A(a;i) • ... • A{xn)- Hence, also 
(III) is satisfied. 
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FIG. 1. Depiction of two-point correlators (ula^) for different point 
separations n in dependence on J/T for antiferromagnetic (ferro- 
magnetic) exchange to tlie rigiit (left) of 0. 



Finite temperature. The coefficients in (j6]l neither depend 
on N nor on the parameters ej. These only enter in lu via the 
dependence of a and G on them. Therefore, taking the Trotter 
hmit N ^ co'is easy. 

The evaluation of the correlation function consists in nu- 
merically solving (jsj for a{x). Then (10 1 is solved for G 
and eventually lo is calculated from (j9i. The definition of 
the coefficients in (j6]) may be found in |4|, but an efficient 
computation is due to 1 12|. Here we use their results for den- 
sity matrices £>„ with n up to 7 sites, and combine this with 
the finite temperature results for the function w. The formu- 
las for Dn{xi, x„) allow for taking the homogeneous limit 
xi = ... = Xn = and yield simple algebraic expressions 
involving uj and its derivatives. 

Just as an illustration of the explicit expressions we like to 
show formulas for the 2-point correlation functions obtained 
from Dn for rt = 2, 3, 4 



KV2^> = ^ (0, 0) , Ka3^) = ^ (0, 0) + ^ (1, 1) - i (2, 0) 



= (0,0) 
(1,0) 
[4(1,1) -2(2,0)] 



^^(M) + ^(2,2) 



27 



(3,1) 



^(1,0) + ^(2,1) + 1(3,0) 



:(3,1) 



1 1 



(2,0) 



1 



+ ^(2,2) (11) 
6 



where {j,k) :— dld20j{xi,X2)\xi=x2=Q- We derived simi- 
lar formulas for the cases n — 5, 8 but these are by far 
too long to be shown here. The correlators {(jfa-n) ana- 
lytic functions of J/T along the entire real axis with zero of 
(n - l)th order at J/T = 0. In Fig.[T] we show the results for 
the antiferromagnetic chain: data at negative values of J/T 
correspond to results for the ferromagnetic chain at \J/T\. In 
the ferromagnetic case the correlations are strictly positive, in 
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3 


4 


5 


6 


7 


0.0470 


0.1070 


0.3268 


0.5014 


0.9013 


1.1957 


1.7761 


1.1283 


0.6419 


0.8714 


0.7521 


0.8652 


0.7971 


0.8699 



TABLE I. The low-temperature expansion coefficients 7^ of the cor- 
relations {(Jiaf^r) for r = 1, 7 (2nd row) and the ratio of 7,. with 
the CFT prediction -y^^'^ = (3rd row). 



the antiferromagnetic case, the correlations are negative (pos- 
itive) for even (odd) n. 

In if TSl the low-temperature behaviour of uj{xi,X2) was 
calculated with the result 



+ (1 - {xi - X2f) cos(7r(a::i + X2)), (12) 



where luq denotes the T = limit of lo in the thermodynamic 



limit. From this expansion and ( 1 1 1 we obtain exphcit results 
for the correlations like 



'1^3/ 



i-^ln2 + 3C(3)- 



TT 

72 



(r/J)2 (13) 



Similar but more lengthy expressions up to ((tJ^(t|) have been 
derived. These results can be written in the following form 



'l"l+r 



)o{l -MT/Jf) 



(14) 



where explicit numbers for are given in table |l] 

For a (primary) field with scaling dimension x the two-point 
coiTelator at distance r is given by 



GriT) = C 



ttT/v 



sinh nrT/v 



2x 



c 



^{■KrT/vY 



(15) 

where the sound velocity of the elementary excitations is 
V — 2ttJ. For the spin-spin correlations we have to set 
X = 1/2. Hence, the CFT prediction of the coefficient in ( 14i 
is 7^^"^ = r^/24. In tablenlresults for the ratio Jr/jr^"^ are 
given. The values are of the order 1, but deviate significantly 
from 1 which we attribute to the simple form of the CFT pre- 



diction ( 14 1 which is strictly valid only for conformally invari- 
ant models without marginally irrelevant perturbations which 
exist in the isotropic Heisenberg chain. Also, it is likely that 
the sequence ^r/lr^^ has two different accumulation points 
for even and odd r, respectively. 

In Fig. [2] we show the entanglement entropies 5„(T) for 
blocks of size n for n = 1, 7. Note that we used the log- 
arithm with base 2 (the total number of local states) in the 
definition of the entropy. Hence the high temperature asymp- 
tote of Sn is identical to the length of the block n. The low- 
temperature limit scales with l/31og2(n). 

Finite Size. Similar to ([T2]| we find in the large L limit 



UJ 



Wo 



+ ^ (1 - (z^i - X2f) cos(7r(a;i + X2)). (16) 
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FIG. 2. Depiction of the entanglement entropies of Dn as functions 
of temperature. 
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FIG. 3. Antiferromagnetic Heisenberg chain at T = 0: Depiction 
of 2-point correlators (afa^) in dependence on chain length L for 
different point separations n. Solid lines are finite size fits. 



For finite size, the physical density matrix is Z?„ = 
Z)„(l/2, 1/2). Hence, all the above formulas for low T 
turn into their counterparts for large L if we replace (2 J/T)^ 
by —L^/tt^ leading to 

(aiVJ+,) ~ KaJ+,)o(l + 47.^Vi')- (17) 

In contrast to finite temperature, finite size increases the cor- 
relations. In Fig. [3] we show data for chains of length L up 
to 128. Note that correlators {afaf^), (crj^cr^) for the same L 
coincide if n + m = L + 2. 

In conclusion we have derived exact results for correlation 
functions of the Heisenberg spin chain for finite temperature T 
(finite system size L). In the conformal regime of low T (large 



L), the corrections to the ground-state results in the thermody- 
namic limit are additive (L^^) terms. The exponents are 
universal and agree with CFT predictions. The coefficients 
for two-point correlators at strictly finite lattice separation 
are non-universal, but correspond precisely to those quantities 
that are of central interest in many, especially numerical ap- 
proaches like the diagonalization of Hamiltonians or Monte- 
Carlo simulations. We also managed to derive multi-spin cor- 
relations and as an example we showed exact data for the 
entanglement entropy for arbitrary temperatures with smooth 
transitions from T = to T = oo with rather different depen- 
dence on the length of the chain segment. The central expres- 
sion for all (static) correlation functions shows a remarkable 
structure, it is a sum of products of nearest-neighbour corre- 
lators. We also explained how the central expression for the 
coiTelations can be derived from a set of 'discrete' functional 
equations for the density matrix. We are convinced that this 
method is very powerful and applicable to other seminal mod- 
els with higher spins or different symmetry groups as well. 
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